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Abstract. We determine rates of convergence of rank-based interacting diffu- 
sions and semimartingale reflecting Brownian motions to equilibrium. Conver- 
gence rate for the total variation metric is derived using Lyapunov functions. 
Sharp fluctuations of additive functionals are obtained using Transportation 
Cost-Information inequalities for Markov processes. We work out various ap- 
plications to the rank-based abstract equity markets used in Stochastic Portfo- 
lio Theory. For example, we produce quantitative bounds, including constants, 
for fluctuations of market weights and occupation times of various ranks for in- 
dividual coordinates. Another important application is the comparison of per- 
formance between symmetric functionally generated portfolios and the market 
portfolio. This produces estimates of probabilities of "beating the market". 



1. INTRODUCTION 

Informally a rank-based model is a multidimensional Markov process whose in- 
stantaneous dynamics is a function of the order in which the coordinates can be 
ranked. When the process is a diffusion its movement can be described by a stochas- 
tic differential equation (SDE) as follows. Let n be the dimension of the process. 
Let Si and cr,-, 1 < i < n, be two finite collections of real and positive real con- 
stants, respectively. For any vector x e E n , let xm < Xf 2 ) < • • • < X( n ) denote its 
ranked coordinates in increasing order. Consider the following system of stochastic 
differential equations: For each i = 1,2, ... ,n we have 

(!) dX &) = (l> ' l {x t [t)=x u) (t)}\ dt + (l> • ^x^x^yj dW W- 

Here W = (Wi, . . . , W n ) is a system of jointly independent one-dimensional stan- 
dard Brownian motions. Assumptions and conditions guaranteeing the existence of 
such processes will be discussed later in the text. 

Different versions of the particle system in ([I]) have been considered in several 
recent articles. Among the more recent ones, see Banner, Fernholz and Karatzas [4], 
Banner and Ghomrasni [5. , McKean and Shepp [30] , Pal and Pitman [32] , Jourdain 
and Malrieu [53], Chatterjee and Pal [HJ[T2], Ichiba and Karatzas [52j, Ichiba et 
al. [53] . Pal and Shkolnikov [53], and Shkolnikov [101 E]. Related discrete time 
processes are studied in the context of competing particle systems by Arguin and 
Aizenman Ruzmaikina and Aizenman [38], Shkolnikov [39], and Racz [35j . We 
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refer the reader to the above articles for the full list of applications of such processes 
(keywords: McKean-Vlasov equations, competing particles, dynamic models of spin 
glasses, models of equity markets, and queueing theory). 

Processes satisfying ([!]) do not have any equilibrium. For example consider the 
average of all the coordinates. The resulting process is a linear diffusion with a 
constant drift and a constant diffusion coefficient. This process clearly has no 
equilibrium. The interesting long-term asymptotics arise when one considers the 
centered process by subtracting the average from each of the coordinates. Under 
appropriate conditions (see below) , the centered coordinate process is positively re- 
current in the neighborhoods of the origin and has a unique stationary distribution. 

Another interesting and useful feature of this model is the process of spacings: 

Yi(t) = X {i+1) (t) - X w (t), i — 1,2,. ..,71 1. 

The vector Y = (Yj., Ya, • ■ ■ , Y n -\) has the law of a semimartingale reflecting Brow- 
nian motion (SRBM) in the positive orthant (R+)"'~ 1 with a constant drift and 
diffusion coefficient and a constant oblique direction of reflection on each face of 
the orthant. This is an important class of processes that arise as a heavy-traffic 
limit of queues. Please see the survey by Williams [13] . 

A particularly interesting example of the rank-based models is the so-called Atlas 
model in which one takes Ui = a for all i and Si = for all i > 1. At any moment, 
the smallest particle gets an upward drift of 8\ = 6 > while the others are 
locally independent Brownian motions. The asymmetry of the interaction among 
the particles sets it apart from the usual models of colliding particles (see, e.g., 
Harris [5U], Swanson [12] )• In discrete time, this model can be compared with the 
asymmetric exclusion process or the Hammersley-Aldous-Diaconis process. 

In this paper we shall determine exponential rates of convergence to the equilib- 
rium for rank-based models and the derived reflecting Brownian motions. We also 
obtain sharp Gaussian fluctuations around the equilibrium mean for additive func- 
tionals of these processes. The main difficulty one encounters with these processes 
is the lack of smoothness in its drift and diffusion parameters. Classical theorems 
often cannot be used and novel methods have to be invented. Our approach is 
a combination of several tools: an extended construction of Lyapunov functions 
for SRBM introduced in Dupuis and Williams [15], recent advances in Transporta- 
tion Cost-Information inequalities as in Bakry et al. [3] and Guillin et al. [115] and 
classical Poincare inequalities satisfied by the associated Dirichlet forms. 

We should add that geometric ergodicity for certain recurrent SRBMs has been 
shown in a recent paper by Budhiraja and Lee [9]. However, these results lack 
constants and the proofs cannot be modified to yield explicit estimates. Our results 
have explicit constants wherever possible. We also give a shorter and simpler proof 
of geometric ergodicity itself. 

One of our primary motivations is to solve certain open problems described in the 
survey by Fernholz and Karatzas [IB] . These are related to the area of Stochastic 
Portfolio Theory to which we provide a very brief introduction. 

1.1. A brief introduction to stochastic portfolio theory. Stochastic Portfolio 
Theory (Fernholz [T7], Fernholz & Karatzas |T5] ) is a descriptive theory of equity 
market models (i.e., a dynamical model of total wealth that various companies 
raise through their stocks) which aims to be compatible with data on long-term 
market structure. This is a departure from the usual Black-Scholes models that 
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FIGURE 1. Capital distribution curves: 1929-1999 



are normative and are not supported by data. One significant difference between 
the two models is that while the Black-Scholes model assumes the principle of 
no-arbitrage in all its formulas, models in SPT, in fact, try to uncover arbitrages. 

The importance of the rank-based models stems from the fact that they match 
the data of the capital distribution curve. To explain this, suppose index i denotes 
a typical company listed in any of the major U.S. stock exchanges. Let Si(t) 
be its market capital (i.e., number of shares x price of each share) at time t. Let 
S(t) = (Si, S2, ■ • ■ , S n )(t) denote the vector of market capitals of all such companies. 
Here n is of the order of several thousands. We define the market weights as 

_ Sj(t) 

A 4 * W — sr^n q~T7\ ' * L, Z, . . . ,Tl. 

It has been observed for over eight decades that if we arrange the market weights in 
decreasing order then they display a power-law decay. In other words, the ith largest 
market weight is proportional to i~ a . See Figurejljwhich plots log ^t( n+1 _j) (t) versus 
log? for one instance of t per decade. These curves are called capital distribution 
curves. It is known (see [I], [TT]) that if we define Xi — log Si to be a rank-based 
model, then, under appropriate assumptions, such a power law decay can be proved 
in equilibrium for large n. 

However, the stability of the market weights over time is more mysterious. It is 
clear that market parameters do not remain constant. In fact, the number of listed 
companies have grown ten fold. The reason for this stability has three reasons: 
(i) the limiting shape is independent of the parameters satisfying certain weak 
conditions, (ii) a fast return to equilibrium when perturbed in the parameters or 
dimension, (iii) tame fluctuations during the period that it takes the process to 
mix. (i) has been shown in [11 . In [33], the authors have proved a fluctuation of 
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order vT for the shape of the market weights over intervals of time of length T. 
This shows (iii). In this article we take up (ii) and show an exponential rate of 
convergence for market weights. 

Another important topic of analysis is the performance of portfolios when the 
market is modeled as rank-based. A portfolio (jri(t), . . . , 7T n (t)) is a random process 
that takes values in the simplex {x S R" : Xi > and Xi = 1}. The quantity 
7Tj (t) represents the proportion of wealth an investor invests in stocks of company 
i at time t. When tt = fj,, the market weights, the portfolio is called the market 
portfolio. The latter is of central importance since the market portfolio determines 
various index funds (i.e., portfolio tracking an index such as S&P 500). An impor- 
tant question in practice is the relative performance of portfolios when compared 
to the market portfolio. 

Let V n (t) denote the wealth of an investor who has followed portfolio tt with 
an initial investment of $1. Often tt = H(/i) is a function of the market weights. 
Such portfolios are called functionally generated (see the survey [IB])- We focus 
on estimating tail probabilities of the quantity V" (t) /V^(t) and its reciprocal for a 
functionally generated portfolio tt. These are probabilities that the portfolio tt will 
beat the market which determine its attractiveness to investors. 

This is closely related to finding relative arbitrage opportunities. According 
to the definition in |18j . the rank-based models are not diverse, i.e., the market 
entropy process is not bounded from below. This implies that the entropy weighted 
portfolio (see Section |4~Tj below) docs not necessarily provide a long-term arbitrage 
opportunity over the market portfolio. Thus a natural question is an estimate of the 
tail probability of the ratio of value processes between the functionally generated 
portfolio and the market portfolio. 

1.2. Main results. Define the set I as {1, ... , n} and consider the SDE in ([lj: 

(2) dXi{t) =J2 5 r l {x l {t)=x (0) (t)} dt + a i ■ l {x l {t)=x U) (t)} dw ^ t )> { e L 

jei jei 

The weak solution of ([2| exists and it is unique in law (e.g., [4], [7]). 

We assume the following constraint on the drift constants: for 1 < k < n — 1, 
the constants 

k 

(3) a.k ■= 2^~^ (#j — S) are strictly positive, 

i=i 

where 5 := J2jei <W n ^ s aver age drift. This is a stability criterion that guar- 
antees the existence of an invariant distribution for spacings or "gaps" (see [2~j]V 
For some of our results we will also assume the condition 

/as 2 2 2 2 

(4) °2 - CT i = ■ ■ ■ = <r n ~ a n-l 

on the diffusion coefficients. This is an assumption that is directly motivated by 
data on volatility. See Figure 13.6 in [TH] which shows estimated volatilities to be 
almost linearly decreasing with rank. 

Let v denote the law on M™ -1 of an (n — 1) -dimensional random vector with 
independent coordinates, where the fc-th coordinate is distributed according to the 
Exponential distribution of parameter 

&k := 2a k (al + crl^y 1 . 
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It is shown in Corollary 2 of [24] that under ^ the distribution of the spacing 
system (Y(t) := (Yi(f), . . . ,Y n _i (*)),* > 0) consisting of 

^(t) := - X (j) (t), j = l,...,n- 1, 

converges in total variation norm as t — > oo to its unique stationary distribution. 
Moreover, under condition Q the latter is given by the measure v defined above, 
as explained in section 5 of |24| . We call a probability measure k supported by 
(R + )™ -1 a spacing distribution. The above v is a spacing distribution. 
In this setting we obtain the following theorem. 

Theorem 1. Let condition ^ hold and suppose that = 1, i £ I. Then for every 
initial spacing distribution k that is absolutely continuous and such that dn/dv is 
in L ((R_|_) n_ , v), and for every bounded function u of the spacings with v(u) :— 
Jm \ n -i u dv = and Var u (u) := J™ u 2 dv — a 2 , and t,r,e > 0, the following 
estimate holds: 
(5) 

t 



P 



u(Y(s))ds > r 



< 



d,K 
dv 

Here S(u) 
(0) 



exp 



t 



13 

sup \u(x) 



\5 2 (u) 

u(y)\ is the range of u and 
4A„ 




mini<Kr 



~2 ' 



It is 
when 



where X n is given in ( 10 ) below. 

Remark 1. Notice the effect of various parameters in the expression (J5j|. 
exponential in t showing a geometric rate of convergence. It is Gaussian in r 
r is large, which shows a strong concentration around the mean. This is evident in 
real data on additive functionals such as cumulative excess growth rate as can be 
seen in Figure 11.3 in [18]. 

Theorem 2. Let conditions ^ and Q be satisfied. 

(i) For every function u £ L 2 ((M + ) n_1 ; v) for which v(u) — 7 ||Vu|| G 
L 2 ((M + )™~ 1 ; v) and the conditions of Lemma ^ below hold, 
variance bound 



obta 



E" 



u(Y(s))ds 



0{t~ 



(ii) Furthermore, for every < e < 2 and every initial spacing distribution k 
that is absolutely continuous with respect to v and such that dn/dv belongs 
to L 2 / e ((M + ) n_1 ; v), the bound of part (i) is modified to 



u(Y(s))ds 



oit- 1 ). 



Remark 2. Under additional assumptions on the function u in Theorem [2] we obtain 
bounds of order i~ p / 2 on the p-th moment of the corresponding additive functional, 
and with more assumptions also the Gaussian concentration property for the latter 
(see Lemma [6] below for the details). Moreover, it is clear from the proof that the 
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result of Theorem [T] applies to any recurrent normally reflected Brownian motion 
in a polyhedral domain with constant drift vector and identity covariance matrix. 
Similarly, the result of Theorem [2] holds for any recurrent semimartingale reflected 
Brownian motion in a polyhedral domain with skew-symmetric data in the sense 
of 03]. 

In the absence of condition Q we view the spacings process as an SRBM in the 
sense of [TS]. The general theory of SRBMs is quite tricky and depends heavily 
of properties of Skorokhod maps. We will use the conditions of Theorem 2.6 in 
|15j which yield the existence of a unique invariant distribution. In technical terms 
these conditions are (i) a completely-S condition on the reflection matrix, and (ii) 
the statement that all solutions of an associated deterministic Skorokhod problem 
are attracted to the origin. We reprove the geometric ergodicity result of Budhiraja 
and Lee [9] with a much shorter and simpler proof. 

Theorem 3. Let Y be an SRBM satisfying the conditions of Theorem 2.6 in [15] . 
Denote by P t (y,.), y € t > the transition probabilities of Y and let 

1 1. 1 1 TV be the total variation norm. Then there exists a measurable function M : 
(R+)™ _1 — > [0, oo) and a constant < £ < 1 such that it holds 

(7) \\Pt{y,-)-v\\TV<M(y)C 

for all y £ (K+) n_1 , t > 0, where v is the unique invariant distribution ofY. 

1.3. Outline. The article is organized as follows. Section [2] establishes rates of 
convergence for SRBMs and their additive functionals under the condition of skew- 
symmetry. In Section [3] we prove exponential rate of convergence for a general 



SRBM under a technical condition on the Skorokhod map. Section 4.1 establishes 
comparison bounds on the performance of functionally generated portfolios with 
respect to the market portfolio. Finally, in Section |4.2| we give estimates on the 
fluctuations of the market weights and of the time spent by a given market weight 
in a given rank. 

2. Convergence to equilibrium under the skew-symmetry condition 

To be ready to prove Theorem [TJ we start with some preliminaries. In the first 
subsection of this section we recall the construction of certain normally reflected 
Brownian motions with constant drift. In the following subsection we explain how 
the latter are related to the spacings process in the case that cr* = 1, i € / and give 
the proof of Theorem [T] in this case. Finally, in the last subsection of this section 
we extend our approach to the setting of Theorem [T] 

2.1. Normally reflected Brownian motion with constant drift. Fix a di- 
mension neN, define the open wedge 

H = |a; € K n : ^.t^ > 0, and x x < x 2 < ■ ■ ■ < £ n | 

and let 7 = (71, ... ,7 n ) be a given vector of constants. 

We consider a normally reflected Brownian motion (RBM) on H with a constant 
drift vector 7, an identity covariance matrix, and normal reflection. This process 
can be obtained by the following Dirichlet form construction as described in section 
2 of Burdzy et al. [TU]. 
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First, define the positive measure 

m(dx) :— 1h(x) exp(2(7, x))dx. 

Next, let T be the collection of continuously differentiable functions on H such that 
both the function and its gradient are square integrable with respect to m. That 
is, 

(8) T = {ueL 2 (H;m) : Vti € L 2 (H; to)} . 

Over this domain we define the symmetric Dirichlet form 



£(u, v) := — I (Vu(i), Vt'(i)) m(dx), for ti, w £ J. 
2 J H 



It can be shown (see [TO]) that (£,.F) is a regular Dirichlet form on h 2 (H;m) in 
the sense that C C (H) n J- is dense both in C C (H) in the uniform norm and in J- 
with respect to the Hilbert norm 



\/£i(u, u) :— ^£{u, u) + (u, u) 



L, 2 (H;m)- 



Hereby, as usual, C C {H) denotes the space of continuous functions with compact 
support in H . 

A strongly continuous Markov process X is called an RBM on H with a constant 
drift 7, if it is symmetric with respect to the measure to and the associated Dirichlet 
form is given by {£ , J 7 ) in the sense 

1 

t->Q t 



T = -j u e L (H\ m) : Jim - {u — P t u, u) L 2( H . m) < oo 



£(u,v) = lim - (u- Ptu,v) L 2 {H . m} , for u,v € T. 

Here (Pt)t>o refers to the transition semigroup of the process. 

The function h{x) = 1h{x) exp(2 (7, x)) is the density of an invariant measure for 
X. When h is integrable with respect to the Lebesgue measure on H, the normalized 
measure (which we will continue to denote by m) constitutes the unique invariant 
distribution of the RBM. 

It is known that this process is unique in law and admits a semimartingale 
decomposition, given by the solution of a deterministic Skorokhod problem applied 
to the path of a Brownian motion with drift vector 7. In particular, the process 
has a local time at the boundary and, started from any point in H, almost surely 
never visits the origin [43 . 

2.2. Proof of Theorem [TJ We now return to the rank-based model described in 
the introduction and assume throughout this subsection that <7j = 1, % € I. As 
explained in the proof of Theorem 8 in [35], the law of a suitably shifted vector 
of ordered processes in this model is identical to that of an RBM in H. More 
specifically, let j3 denote an independent one-dimensional RBM with a negative 
drift -9 (9 > 0). Define 

Yi(t) := X (t) (t) - X(t) + i = 1,2, . . . ,n. 

Note that £ i6/ ^(t) = \fnl3{t) > for t > 0. In addition, for every t > the 
spacings between the components of the vector Y(t) are the same as the original 
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spacings, and, moreover, the process Y is an RBM in H, as described above, with 
a drift vector 7 determined by 

n , _ Q \ l n ~ l Q n 

(7, y) = [ Si ~ 5 ~ ~r ) & = ~ o ak ^ k + i ~ y fe ) ~ 7= y e H - 

It follows from |32|, Lemma 9] that the corresponding invariant measure m is in- 
tegrable with respect to the Lebesgue measure on H and that, under its normalized 
version, the law of the spacings {jjk+i — Vk, k = 1, . . . , n — 1) is given by a product 
of Exponential distributions with parameters ak, k = 1, . . . , n — 1. 

For any y S H , we make the transformation 

yk = Vk+i-yk, fc = 1, 

which maps configurations of ordered particles to the corresponding vectors of spac- 
ings. Next, we let T be the unique linear extension of this map to the whole of 
R". The push-forward of the measure m by T, which is a probability measure on 
the positive orthant will be denoted by v and is given by the product of 

Exponential distributions with parameters a^, k = 1, . . . ,n — 1. 

Consider any function / : (]R + )™ _1 — > R which is continuously differentiable. The 
gradients V y / on and V^(/ oT) on H satisfy the simple linear equation 

d(f o T) _ ^ Of dy, _ ^ a/ 



<9yj ^ dyj dyi 13 dyj ' 

where A = (ay) is an n x (71 — 1) matrix given by 

if l<j = i<n — I, 
if 1 < j = z — I < n — 1, 
otherwise. 

Moreover, it is clear that / and its gradient are in L 2 ((R+)™ -1 ; v) if and only if 
g := f o T and the gradient of g are in li 2 (H; m). In fact, we have the following 
bound: 

(9) f IIVJII 2 dv < A„ / ||V s3 || 2 dm. 

i(R + )™~ 1 

Here, 1/A„ is the minimal eigenvalue of the matrix A' A, with A' denoting the 
transpose of A. To compute the latter eigenvalue explicitly, we note that 2- Id— A' A 
is the (n — 1) X (n — 1) tridiagonal matrix which has zeros on the diagonal and ones 
next to the diagonal, with Id being the (n — 1) x (n — 1) identity matrix. It is well- 
known (see e.g. [5S] and the references there) that the eigenvalues of a tridiagonal 
matrix of this type are given by 2 cos— , k — l,...,n — 1, so that the largest 
eigenvalue of 2 • Id — A' A is given by 2 cos - . Thus, 

^ A «=2-L- 

n 

In particular, we see that X n grows quadratically in n. 

Now, recall that a probability measure /j, on R" is said to satisfy the Poincare 
inequality with the Poincare constant Cp > if for every continuously differentiable 
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function / : R n -> K. it holds 

(11) Var M (/)=/ fdn-([ fdS <C P ( HV/H 2 ^, 

JR™ Vl» 7 JR™ 

where Var M (/) stands for the variance of / under \i. 

We proceed by recalling the fact that the Exponential distribution with param- 
eter 1 satisfies the Poincare inequality with a constant Cp = 4. This well-known 
result can be found in several sources, such as the book by Ledoux (Lemma 5.1 of 
|29j ) or the article by Barthe and Wolff [6], who also prove several other interesting 
results about the Gamma family. 

It follows by a simple scaling argument that the Exponential distribution with a 
parameter A > satisfies the Poincare inequality with the Poincare constant Cp = 
4/ A 2 . In addition, we use the fact that the Poincare inequality has the following 
tensorization property: if fix , . . . , fj, n are probability measures on Euclidean spaces 
satisfying the Poincare inequality with the constants C\ , . . . , C n , then the product 
measure fj,± X . . . X fi n satisfies the Poincare inequality on the corresponding product 
space with the Poincare constant maxi = i r ...„ Cj. 

The above considerations lead to the following lemma. 

Lemma 4. The measure v defined in the beginning of this subsection satisfies the 
Poincare inequality with the Poincare constant 



C v := 4 • max a h 2 — 4 • ( min ai ) 

Kk<n-1 K \l<fc<n-l / 



Combining Lemma [4] with the comparison bound in ^ and the equation (10), 
we obtain for all functions / S T the estimate 

(12) Var m (/) < 1 • - , • £{f, f) =: 0£(f, f). 

2 - 2 cos - mmi< fe <„_i a% 

Here, the pair (£,J-) is defined for the RBM Y according to the previous subsection. 

We can now finish the proof of Theorem [T] in the case that <7j = 1, i € I. Indeed, 
applying Theorem 3.1 from Guillin et al. 19J, we obtain the estimate of Theorem 
[l]for any bounded measurable function of the particle configuration given by the 
components of the process Y. In particular, the statement of Theorem [l] holds for 
any bounded measurable function of the vector of spacings, as claimed. 

2.3. Proof of Theorem [2} We now give the proof of Theorem [2] In this part 
we assume the condition Q on the drift coefficients and the condition Q on the 
diffusion coefficients. The proof is broken down in several steps. 

Step 1: Reduction to SRBM. The ordered particle system corresponds to the 
following SDE: 

(13) dX(i)(t) = Sidt + <j t dBi{t) + ^dLi- 1:i (t) - ^dL i:i+1 (t), i = l,...,n, 

where refers to the local time of collisions between the i-th and the (i + l)-st 

ranked particles (see section 3 of [!])• In other words, this is the semimartingale 
local time at zero for the process Xu + i)(t) — Xu^t), t>0 normalized according to 
the Ito-Tanaka formula. 

In the following we closely follow the analysis and linear algebra done in Ichiba 
and Karatzas [531 Section 3.2.1]. The process (Xn), . . . , X^) is a semimartingale 
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reflected Brownian motion (SRBM) in the sense of [43], with a constant drift vector, 
a constant diffusion matrix and normal reflection, taking values in the Euclidean 
closure of the wedge 

{(xi,...,x n ) G E™ : Xx < ... < x n }. 

Next, we consider the process of spacings Y, whose components are governed by 
the SDEs 

dYi(t) = (5 i+ i - 5i)dt + a i+l dB i+ i{t) - aidB^t) 

+dL i}i+1 (t) - ~dL i+14+2 (t) - ^dL l _ 14 (t), 

i = 1, . . . , n — 1. It is easy to verify that Y is also an SRBM in the sense of [33], 
taking values in the positive orthant (R+)"™ 1 , and having a constant drift vector 
7, a constant diffusion matrix S and a reflection matrix R. 

More specifically, the drift vector is given by 7 = (62 — Si, ... , S n — <5„-i). The 
covariance matrix S = (£ij)i<i t j< n — 1 is tridiagonal with components 



(14) 



'of + of +1 , when l<i=j<n — 1, 

-erf, when 2<i = j + l<n, 

-of 1 15 when 1 < i = j — 1 < n - 

0, otherwise. 



The reflection matrix R = (rij)i<ij< n -i is also tridiagonal and its components are 
given by 



(15) 



1, when l<i=j<n— 1, 

1 when 2<i=j + l<n, 

when 1 < i = j — I < n — 1, 
0, otherwise. 



2 ' 
1 

2 • 



Let E denote a nondegenerate square root of the matrix S, that is, an invertible 
(n — 1) X (n — 1) matrix which satisfies EE' = S. Note hereby that the matrix S is 
stricly positive definite, since none of the processes <Ji + iBi(t) — Oj-Bj(t), t > for 
i = 1, . . . , n — 1 can be obtained as a linear combination of the other n — 2 processes 
of the same type. 

Now, one can make the transformation Z{t) := E _1 y(t), t > to obtain an 
SRBM in a polyhedral domain G given by the image of the orthant (]R + ) ,l_1 under 
E _ . This new SRBM has a constant drift vector given by E _1 7, identity covariance 
matrix, and a new reflection matrix *H := E _1 _R. 

Hence, for the sake of the proof of Theorem [2] we can restrict our attention to 
the case of an SRBM taking values in the closure of a polyhedral domain G with 
data of the following form: 

(i) drift vector E _1 7 and reflection matrix 9^, 

(ii) a cone G of dimension n — 1 given by the dual description: 

G = {x: (Ex), > 0, i = 1, . . . , n - 1} . 

In particular, the (n — 1) x (n — 1) matrix N whose columns are the unit 
normal vectors on the faces of G is given by the normalized row vectors 
of E. In addition, the equality S = EE' implies that the norm of the i-th 
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row of £ is equal to Thus, if we let D denote the diagonal matrix 

comprised of the diagonal elements of 3, we have N = Yl D~ x l 2 . 
(iii) Due to the condition Q the reflection matrix £H admits the decomposition 

D\ = N + Q, 

where the diagonal elements of N'Q are all zero. In fact, under Q the 
so-called skew- symmetry condition 

(16) N'Q + Q'N = 

holds (see page 25 in [23] for details). 

Following Williams [33] we call the process Z an SRBM associated with the 
parameters (N, Q, S _1 7). Such an SRBM has the following semimartingale decom- 
position: 

(17) Z{t) = Z(fl) + ITS'* + W(t) + (N + Q) L(t), 

where L is the vector of the accumulated local times at the faces of the boundary 
of G. 



Step 2: Application of duality. Next, we define 7 := 2(7— N~ 1 Q) S _1 7 
and let 7 denote the vector 7 — S _1 7. The above considerations show that the 
following result from Williams [J3J Thm 1.2, Cor 1.1] applies in our setting. 

Theorem 5. Consider the measure p on G whose density with respect to the 
Lebesgue measure on G is given by 

p(x) = exp((7,x)). 

The SRBMs associated with (N, Q,X _1 7) and (iV, — Q,j) are in duality relative to 
p and p is an invariant measure for these two processes. In particular, if p is finite, 
the normalized measure (which we will also denote by p) is the unique stationary 
distribution for each of the two processes. 

The duality being referred to is in the following sense: let {Pt)t>o and (Pt)t>o 
denote the transition semigroups for the two SRBMs in the theorem. Then for all 
continuous functions /, g with compact support in the closure of G we have 

(18) / (P t f)(x)g{x)p(x) dx = f f(x)(P t9 )(x)p(x) dx. 

JG JG 

Now, let C^°(G) be the collection of infinitely diffcrcntiablc functions with com- 
pact support in the open cone G and let C and C denote the respective gener- 
ators of the two SRBMs in the statement of the latter theorem on the domain 
C C °°(G) C L 2 (G;p). That is, 

(19) £=(E- 1 7 ,V) + iA, £=(7,V) + ^A. 
From the said duality we conclude: 

(-Af,<7> L2 ( G;p) = (-4,/) L , f,geC™(G). 
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Next, for f,g £ C£°(G), we define the symmetrized Dirichlet form 8 a with 
domain C™{G) by 

(20) 



1 



1 



2 \ / L 2 (G;p) 



Moreover, using integration by parts we get 

/-Cg-Cg,f) = f h( 7 ,(V ff )(x))-(A 3 )(x)]/( a ;)exp(( 7 ,x))dx 

' L 2 (G;p) Jq 



(21) ^ Jg 



G 



Thus, 

(22) £°(f,g) = l -f ((Vg)(x),(Vf)(x))exp({j,x)) dx. 

•J G 

In other words, £" is the pre-Dirichlet form for an SRBM with a constant drift 
vector, identity covariance matrix and normal reflection on the faces of G. It is 
known in such case (see the references in Burdzy et al. [TO]) that the Dirichlet form 
is closable in L 2 (G; p) and that the resulting closed extension is regular in a sense 
made precise in the previous subsection. Moreover, the domain of the Dirichlet 
form E a is given by ([8]), with H replaced by G and m replaced by p. 

Step 3: Forward-backward martingale decomposition. Now to analyze ad- 
ditive functionals, note that the duality ( 18 ) between (Pt)t>o and (Pt)t>o can be ex- 
pressed in the following way. Suppose {Z(u), < u < t} is an SRBM (N, Q, S _1 7 ) 
such that Z(0) ~ p (which is assumed to be a Probability measure). Define the 
time-reversed process Z(s) = Z(t — s), < s < t. Then {Z(u), < u < t} is an 
SRBM (N, —Q, 7) with Z(0) ~ p. This follows since the time-reversed process is 
obviously Markov, and the duality equation ( 18 ) determines its transition kernel. 
Now let v : G — » K be a continuous function such that 

(23) IMloo < 00 > ano - / v(x)p(x)dx = 0. 



Suppose U : G — > K is a twice continuously differentiable function (i.e., U £ C 2 {G)) 
such that 

(24) - (l + C^j U(x) = 2v(x), x £ G, and (VU(x))' n(x) = 0,i£ dG. 

Here n(x) is the inward normal at the boundary point x. This vector is also one of 
the columns of the matrix N. 

In other words, U is the solution of the Neumann problem: 

(25) - (7, VU(x)) - MJ(x) = 2v(x), x £ G, and (VU{x))' n{x) = 0, x £ dG. 
The solution to this Poisson equation exists due to Theorem 4.16 in [5]. 

Now, fix a t > 0. We apply a forward-backward martingale decomposition. By 
Ito's rule applied to the semimartingale we get 
(26) 

M(i) :=U{Z{t))-U(Z(0))- [ CU(Z(s))ds- I (VU (Z(s)))' (N + Q) dL(s), 
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is the final element of the martingale VU(Z) ■ W in time [0, t]. 

When we reverse time, applying Ito's rule to Z we get 
(27) 

M(t) :=-U(Z(t)) + U(Z(0))- [ £U(Z(s))ds- [ (VU(Z(s)))' (N - Q) dL(s) 

Jo Jo 

is the final element of another martingale. 



(28) 



Adding Q and (|27j) and using (|24j) we get 

M(t) + M{t) 



= / v(Z(s))ds. 



Thus, for any convex function <f> we have 



(29) 



E<t> 



v{Z(s))ds) < - 



l - [Ecj>{M{t)) + E(j> (M{tj) 



Connection to the Poisson equation. 

Lemma 6. Consider the Poisson equation with Neumann boundary condition in 
(|5l. 

(i) I/||V*7|| p := j G || VU{x)\\ p p(x)dx < oo for some p > I, then 



E 



\ j\(Z(s))dsJ <C p t^ 2 \\VU\\ p . 



where C p is some universal constant. 



(ii) // llVt/ll^ := sup^gg ||V£/(a;)|| < oo, then, for any A e R, we get 
'A f* „ , \ / A 2 ||VC/||L\ 



Eexp J v(Z(s))ds^j < exp 



2t 



Hence, for any r > 0, we have 

p(\ [ v(Z(s))ds > r) < exp f I ■ 

Wo -J- \ mult) 

Proof. We start by estimating (M) (t) and We will only consider (M) 

since the other case is symmetric. By Ito's rule, we get 

(M)(t)= [ \\VU(Z(s))fds. 
Jo 

Moreover, for any k > 1, Jensen's inequality implies 

E (M) k (t) = t k E Q J* \\VU(Z(s))\\ 2 ds^ < t k E Q J* || VU(Z(s))\\ 2k ds^j 



1 



= t K - I E\\VU(Z(s))\\ 2k ds = t k \\VU(x)\\ 2k p(x)dx = t k \\VU\\ 2k . 
* Jo Jg 

The second last equality is due to the fact that the process is running in equilibrium. 

Suppose now ||VE/|| < oo for some p > 1. We invoke Burkholder-Davis-Gundy 

inequality (see [27l p. 166]) for continuous local martingales to infer 

E \M(t)\ p < C P E (M(t)) p/2 < C p f/ 2 \\VU\\ p . 
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The same bound for E M(t) proves (i) by (29 1. 

Now assume that |j Vf/|| < oo. Then, by Ito's rule, it follows that almost surely 
(M(t)) < ||V£/|| t. Hence, by an exponential martingale inequality we get that 
for any A 6 R, we have 



Eexp(XM(t)) < exp 



a 2 iiwii^ 



Again obtaining the same bound for M and combining with (281, we obtain the 
first part of (ii). The rest follows by a standard application of Markov's inequality 
and optimizing over A. □ 

Poincare inequality for the symmetric process. For a large class of functions 
v, one can show that the corresponding ||Vf7|| is in L 2 (G; p). To this end, we prove 
first the Poincare inequality in the setting of Theorem [2] 

Lemma 7. Let the conditions ^ and Q hold true. Then for any f in the domain 
of £ a in (22) which satisfies J G f{x)p{x) dx = 0, we have 



f 2 (x)exp((j,x))dx<C P / \\Wff(x)exp((j,x))dx. 

G JG 

Hereby, Cp is given by min < 4 ^ r " — -a* with A„ being the largest eigenvalue o/S" 1 
and Sfe, fe = 1, . . . ,n — 1 being defined as in the introduction. 

Proof. As remarked in the introduction, condition ^ implies that the process of 
spacings Y possesses a unique invariant distribution, so that the density function p 
can be normalized to a probability density function, which we will refer to as p\. 

We now fix a function / as in the statement of the theorem, define a function 
g : (K + )™ _1 — > K by g{y) = and denote the invariant distribution of the 

process of spacings Y by v as before. From the definitions of v and p\ we see 

(30) f f 2 (x) Pl (x) dx = / g 2 (y) dv{y). 

JG J(R + )™~ 1 

By Theorem 2 in |24j the probability measure v is a product of Exponential distri- 
butions with parameters <Sfc, k = 1, . . . , n — 1, so that as in the derivation of Lemma 
[4] we conclude that v satisfies the Poincare inequality with the Poincare constant 
In particular, it holds 

(31) / 9 2 (y) d V {y) < ^ f \\Vg\\ 2 (y) d V {y). 

In addition, employing the chain rule as in the previous subsection (see ([£])), we 
can bound ||V3|| 2 (y) from above by A„ ||V/|| 2 (ST 1 ?/) for all y e (M+)™ -1 , where 



A n is the largest eigenvalue of (££') —S . Putting together (30 1, (31) and the 



latter observation, we obtain Lemma 7. □ 

A variance bound for additive functionals. It is well-known (see Proposition 
2.1 in the notes by Bakry [2]) that a Poincare inequality as proved in Lemma [7] 
implies that for any function / G L 2 (G; p) such that J G f(x)p(x)dx = one has 

(32) / {P t f{x)f p(x)dx<e- 2t ' c - f f 2 (x)p(x)dx. 

JG JG 
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Here (P*)t>o is the Markov semigroup associated with the Dirichlct form £ a in (22 1. 
We show below that the Poincare inequality entails a variance bound for a large 
class of additive functionals via Theorem [2j 

Lemma 8. For any v £ L 2 (G;p) such that J G v(x)p(x)dx — 0, suppose that the 
function 

/>oo 

(33) U{x) = / P t v(x)dt, x e G, 

Jo 

is in G 2 (G) and satisfies the Neumann boundary condition. Then U is a solution of 
the Poisson equation p5| ). Moreover, if ||Vu|| € L 2 (G;/)), then ||VZ7|| G L 2 (G;p). 

Proof. We first claim that [/ is finite everywhere on G. To see this choose < 6 < 
2/Cp and note that by the Cauchy-Schwarz inequality 

2 



(34) 



P t v(x)dt 



< 



e~ st dt 



„<5t 



(P t v(x)) 2 dt 



The quantity on the left is infinite only if the quantity on the right is. But the 
right-hand side is integrable with respect to p{x)dx by (32) and our choice of 8. 

Let A denote the generator of the semigroup (Pt)t>o- Then, by the usual semi- 
group calculus it follows that 

AU(x) = J^ AP t v(x)dt = J (~P t v(x)\dt = -v(x). 

The generator A is an extension of the differential operator (C + £)/2 over G 2 - 
functions satisfying the Neumann boundary condition. It follows that U is a solution 
of the Poisson problem with the Neumann boundary condition. 

For any two points x, y G G, let X X ,X V be two processes with transition ker- 
nels (Pt)t>o driven by the same Brownian motion B while starting from x and y, 
respectively. By an application of Ito's rule 

\Xf -Xy\ = 2 f{X* ~ XI n(Xf))df t + 2 / V, a - ^>ft»)>< + \x-y\, 
Jo Jo 

where £ x and l v are the local time components of X x and X v at dG, respectively. 

Since G is convex, the inner products in the integrands are nonpositive, and hence 

\Xf — Xf\ < \x — y\ with probability one (Lemma 3.1 of [15]). Taking y — > x, this 

shows ||VPtv|| < Pt ||Vt!|| for alH > (see inequality (3.4) in [15]). The rest is an 

application of the Cauchy-Schwarz inequality as in (34). □ 

We can now put together our findings to give a proof of Theorem [2] 

Proof of Theorem 2. Part (i) is a consequence of a combination of Lemmata [6] and 
[8] along with the transformation back from Z toY . Part (ii) is obtained by putting 
together Holder's inequality 



and part (i). 



u{Y(s))ds 



E 1 - 



dn 
dv 



u(Y(s))ds 



2-s-, 



< E 



'(f) 2 V[(i />(.))*) 



r/l 

7 



□ 
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Remark 3. It must be known that U as defined in Lemma [8] is automatically in 
C 2 (G) and satisfies the Neumann boundary condition. However, at the present 
moment we cannot find a suitable reference. 

3. Convergence to equilibrium for a general SRBM 

Throughout this section we consider a general semimartingale reflecting Brow- 
nian motion Y in the (n — l)-dimensional orthant 6 := (M + ) n_1 in the sense of 
(15j . The set of parameters for such a process is given by a drift vector 7 G ]R n_1 , 
a covariance matrix H = € ]R(" _1 ) X (™ _1 ) and a reflection vector field r on the 
boundary d& of the orthant & which is constant along each face. Intuitively, the 
process Y behaves as a Brownian motion with drift vector 7 and covariance matrix 
S in the interior of ©, and is reflected in the direction determined by the reflection 
vector field r, whenever it hits the boundary d&. For a precise definition of the 
process Y we refer to section 2 of [15] ■ In particular, the process of spacings 

(X( 2 ) - . . . ,X( n ) - X(n-1)) 



between the particles in the particle system ( 13 ) is an SRBM in & in the sense of 
[T5] (see section 5 in [21] for more details). 
We now give the proof of Theorem [3] 

3.1. Proof of Theorem [3] Step 1. We start by introducing the differential 
operator 

1,3=1 J 

We will construct a function V : ©\{0} — > [l,oo) which satisfies 
(36) CV < - Cl V + c 2 l Bd 

for some positive constants c\ , ci , d, whereby Bd stands for the intersection of the 
ball of radius d around the origin with the orthant ©. This will be done in step 2. 



In step 3 we will combine ideas from section 5 in [14 with (36) to finish the proof 
of the theorem. 



Step 2. In order to provide a function satisfying (36) we start with the Lya- 



punov function W constructed in |15) . The properties of W which will use are the 



following: 








(37) 


w e c 2 (6\{0}), 




(38) 


W(y) -> 00, \y\ 


— > 00, 




(39) 




->o, |y| 


— > 00, 


(40) 


3c>0: (7,(VW)(y)) < 


-c, y € 6\{0}, and 


(41) 


(r(y),(W)(y)> 


< -c, y G 5©\{0}, 



Moreover is a homogeneous function of degree one. 

Since W is homogeneous, all partial derivatives of the form and Q y ^ y . are 
homogeneous functions of degree zero and negative one, respectively. This creates 
the problem that the second partial derivatives have a singularity at the origin. To 
get rid of this problem we will deform the function W near the origin appropriately. 
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Let <p '■ (0, oo ) — > (0, oo) be a twice continuously differentiable function such that 

(i) the following asymptotic properties are satisfied near infinity: 

(42) lim <j>'{x) = 1, lim (j>"{x) = 0, 

X— ¥OD X— ¥OD 

(ii) and the following asymptotic properties are satisfied near zero: 

(43) lim = lim <jj'(x) = 1. 

x-^0+ X x->0+ 

(iii) Moreover, the global condition <fi'(x) > 0, x €E (0, oo) is satisfied and the 
function d>' is bounded. 



Now, define the function U(y) = <j)(W(y)) on ©\{0} and note the following: 
(VU)(y) = cf>' (W(y))(VW)(y), 

d 2 U „ dW dW , d 2 W 



Since W is positively homogeneous and (|42j) holds, there is a d > large enough 
such that if y ^ B^, we have 



n-l 



(44) 



<j>' (W(y)) > 1/2, and £ & 



(y) <c/4, 



where c is the constant in (40 1 and (41). Hereby, the second inequality is a conse- 
quence of ( 42 1 , the fact that all first partial derivatives of W are homogeneous of 

conjunction i 

(7,Vl%)) < 



degree zero and (39). In conjunction with (40), this also implies that 

c 



(45) 



for all y i B d 



Now, since VW is homogeneous of degree zero, and 0' > throughout, we also 
obtain 

(46) (r(y), VU(y)) < 0, for all y <= 96\{0}. 

Finally, each second partial derivative £ .W . is homogeneous of degree negative 



one. Hence, when y approaches zero, we have from (43) 
d 2 U 



(47) 



dy t dy : 



<V) =0(1) + 



dyidyj 

ha 

4>'(W{y)) 



W(y) 



W(y) 



d 2 W 



dyidy^ 



(y) =o(i). 



The last equality follows from ( 43 ) and the fact that the function 

d 2 W 



W{y) 



dyidy. 



(y) 



is homogeneous of degree zero. 



Now, we define V(y) = e xu ^ v \ y € 6\{0} with a positive constant A to be 
chosen later. An elementary computation shows 

n-l 



f) 2 n \ 2 



' dyidyj 



6\{0}. 
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From p4[), the fact that VU is globally bounded and (45 1, it follows that choosing 



A > small enough, we obtain the following estimates on the complement of in 
e\{0}: 



cv 



A( 7 ,VC/> + ^ 



d 2 U 
dyidyj 



A(VC/)'S(VC7) 



5> 

i,i=l 

> / Ac Ac 

< A( 7 ,VC/) + T < -—. 

On the other hand, since U, its gradient, and its Hessian can bounded uniformly 
on Bd, one can choose a large enough C > such that 

£V < -^V + C, for all y G B d . 

Hence, we see that there exist positive constants c± , C2 such that 
(48) CV <- Cl V + c 2 l Bd 

on 6\{0}. 



Step 3. To be able to proceed as in section 5 of [T3], we first show that Y is 
irreducible with respect to the Lebesgue measure on © and aperiodic, in the sense 
of the definitions in section 3 of [14] . Indeed, it is well-known (see e.g. [36] ) that for 
any fixed time t > the probability that Y(t) will be in the interior of the orthant © 



PCX 
Jl 



{V( s )es} 



ds 



is equal to one. Thus, for any open ball B c 6, the quantity _E 

is positive, since a Brownian motion started in the interior of © has a positive 
probability of hitting B before hitting d&. Moreover, the same argument and the 
Feller property of Y (see e.g. [J3]) imply that Y is aperiodic with respect to the 
set B\ and that B\ is a small set for Y . 

From the latter observations we conclude that it suffices to check the condition 
in part (b) of Theorem 5.2 in p3] to obtain our Theorem [3] (note that we cannot use 
part (c) of the same theorem directly, since we do not know that V belongs to the 
domain of the generator of Y and that it is a Lyapunov function with respect to the 
latter). To this end, it suffices to show that there is a function £ : (0, 1] — > [0, oo) 
with £(1) < 1, a petite set B of the process Y and a constant C3 > such that 



(49) 



V(z)Pt(y,dz)<mV(y) + c 3 l B (y) 



holds for all t € (0, 1], ye ©\{0}. 

To show (49), we take the expectation in the change of variables formula in the 



form of equation (3.1) in [U], and use the inequality (46) to obtain 



(50) 



E 



e*V(y(t)) <V{y) + E / e s {CV){Y{s)) ds 



whereby we start the process Y in y. Note that the integrability of the terms inside 
the latter two expectations follows from the fact that the functions V and CV are 
bounded on compact sets and grow at most at an exponential rate as \y\ tends to 
infinity. 
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To finish the proof, we plug in (48) into (501 and use Fubini's Theorem to 
conclude 



(51) 



V(z)P t (y, dz) < e-V(y) + c 2 / e s -*P s (y, B R ) ds 



It remains to observe that our inequality (51 ) is precisely the inequality (33) of [14] . 
so that we can obtain the desired inequality ([49]) from the inequality (51 1 in exactly 



the same way as part (d) of Theorem 5.1 infill is obtained from the inequality 
(33) there. □ 



We now attempt to explicitly describe a Lyapunov function for certain rank- 
based models. Consider the process of spacings (X(2) — -^(l)i • ■ • > ~ X(n—i)) 



in the particle system (13). It is known (see Corollary 2 in [24]) that under the 
conditions 



(52) 
(53) 



Si > 5 2 > ■ ■ ■ > S n , 

Si > Si + i for some i G {1, . . . , n} 



the process of spacings is a positive recurrent SRBM. This includes, in particular, 
the so-called Atlas model introduced in |17j . in which Si > Q, 62 = ■ ■ ■ = 8 n = 0, 
<j\ = . . . = a n = 1. Under the above conditions we can provide an explicit function 
V, which may serve as a Lyapunov function in the proof of Theorem [3] Indeed, 
picking a vector v £ & and setting Vo(y) — e^ v,y ', y € & one has Vq > 1 and 



_^ n — 1 n—1 

ij = l i=l 



This shows that, if one can find a. v E & with (7, v) < and (r(y), v) < 0, y € 96, 



then for a c > small enough the function V(y) = e 



c(v,y) 



y e 6 will satisfy 



(54) 



(CV)(y) < - Cl V(y), y G 6, (r(y), (W)(y)> < 0, y G 8& 



for some constant ci > and may serve as a Lyapunov function in the proof of 
Theorem [3) 

In general, the existence of such a vector v can be easily verified via the Farkas's 
Lemma (see e.g. [371 page 200]) which states the following. There is a vector v with 
(7, v) < and (r(y), i?) < 0, y G 96 if and only if there is no solution to Rx + j = 
in 6. Hereby, i? is the matrix whose columns are the vectors of reflection, which in 



our case is given by ( 15 ). Since R is of full rank, we infer that there is a v satisfying 
our requirements for a Lyapunov function if and only if — i? _1 7 ^ 6. 

Under our condition ( 52 ) one verifies for n > 4 that the vector v — (v\ , . 
with components 



1 



+ e, i — 1, . . . , n — 1 



has the desired properties for any < e < — — — 2^ by using the explicit 
formula for the reflection matrix R given in ( 15 ) and the concavity of the function 

x 1 — y — — . For many other examples of linear Lyapunov functions see |13j . 
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4. Applications 

4.1. Portfolio performance. For the rest of the paper we set = 1, i £ I 
without further notice. As an application of Theorem 1 we shall consider an abstract 
equity market model (Si(-), i £ I) where the market capitalization Si(t) of company 
i £ / at time t > is given by Si(t) :— exp(Xj(t)). That is, each Xi(t) in |2]) gives 
the logarithmic capitalization of company i £ I at time t > 0. This market model 
was introduced by Fernholz in the book [17] and investigated further in the articles 
@], Q2], E2, Q3] and gl] among others. 
By Ito's formula and ([2| we have 

dSi{t) = S t (t)(Y,Sj ■ l{X t (t) = X (j) (t)}dt + dWi(t)) + ^Si(t)dt, 

= £ • = + ^) d< + £ • 

Another application of Ito's formula shows that the market capitalization weights 
fii(t) := Si(t)/Y, je i * e J satisfy 

d/j,i(t) = ( finite variation part ) + /Ltj(t) ((5^ — /ij(i))dWj(t), i £ I , 

where <5,-j is the Kronecker delta. Hence, the corresponding cross variation processes 
grow at the rates 

(55) ^^(i) = Mi(«)w(t)E(^-^W)(^-^W)' M e/2 ' 

fee/ 

For notational simplicity we will from now on write Di for the partial derivative 
with respect to the i-th variable and Dij for the second partial derivative with 
respect to the i-th and the j-th variables. 

A portfolio (7Tj(-), i £ I) in the market above is defined as a stochastic process 
adapted to the Brownian filtration and such that 

7Tl(t) + ---+7r n (t) = l, <>0. 

Its value process V^(-) is defined as a solution to 

For example, 7Tj(.) = Mt(-), i6l gives the market portfolio and we will write V"^-) 
for its value process. 

Following Fernholz |17j , we introduce the family C of portfolio generating func- 
tions G : U — > (0, oo ), which are functions of class C 2 on some open neighborhood 
U C R n of the simplex 

A" := {x = (xi,..., x n ) £ [0, 1]" :x 1 +-»+x n = i} 

and such that the mapping x — > :Ejk>gG(x) is bounded on A™ for all ?' € /. 
According to Theorem 3.1.5 of [T7], given a portfolio generating function G £ C. 
the functionally generated portfolio (7Tj(.), i£ I), defined by 

7r ( (t)= (AlogG(Ai(t)) + l-5^A*j-DjlogG(^(*)))A**(*). * > 0, * G I 
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has a value process V 71 '(•) which satisfies the master formula 

(56) log(V*(t)/V»(t)) = log(G(/i(t))/G( M (0))) + / Q (s)ds, t > 0, 

Jo 

where the drift part is 

(57) .W-iQ^ £ 



Substituting (55) into this formula, we obtain g(-) as a function of the market 
weights /Ltj(-), i € J: 

s(-) - 2G 7 1 n , E AiG(/x(-))^(-)^(-)E(^-^(-))(^-^(-))- 

(j,i)e/ 2 fee/ 

The relative value process V" (•) /V*(-) is determined by the drift process 0. 

We shall now rewrite the drift process g in terms of the spacings process Y of 
Theoremlll The ranked market weights fin_-\(t) < (i(2\(t) < • • • < (ir n \(t) at a given 
time t >0 are obtained from the corresponding value of the spacings process by 
the formula 

M(i)(*) = M i (F(t)), j e/, 

where the functions 1VL, : [0, oo) n_1 — > (0, 1) are defined by 



Mi(yi,...,y n _i) := [1 + 



e »i +e »i+!/2 + ... + e »i+-»»-ii-i. 



(58) M fc (y 1) ... ) y n _ 1 ) : =e«+-+^M 1 ( !/1) ... 1 | /B _ 1 ), fc = 2,...,n. 

For notational simplicity we will write M(y(i)) or M(t) for the value of the process 
(/!(!)(•), . . . , /i( n )(-)) & t a given time t > 0. Let us also introduce the partitions 



iei lei 



such that for every (i, Z) S J 2 and every x = {x\, . . . , x n ) € the component Xi 
is the Z-th ranked from the bottom in the set {x%, . . . , x n }. 

We recall that the domain of the portfolio generating function G is given by A™ 
rather than the set 

(59) A^< := {(xi, . . . , x n ) G A" : Xl < x 2 < ■ ■ ■ < x n }. 
However, if the function G is permutational invariant, in the sense that 

(60) G(xi, . . . ,x n ) = G(x p(1) , . ..x p(n) ) 

for every permutation (p(l) , . . . , p(n)) of /, then we may and will view G as a 
function on A™ Next, let us introduce the functions 

(61) 9jk(x)-= D M G{x)'l Q w{x)-l Q m(x), (j,k)el 2 

(h,i)ei 2 

defined on the simplex A". 

Important examples of portfolio generating functions are: 

• G(.t) = (J2iei x i) 1 ^ p f° r some < p < 1 (diversity) or 

• G(x) = 1 — g J2iei( x i ~ n^ 1 ) 2 (quadratic Gini coefficient) or 

• G(x) = (1 — p)^ 1 log(^ig/ xf) f° r some p ^ 1 (Renyi entropy) or 

• G ( x ) = -J2 t ei x ^°S x i (entropy) or 
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• G(x) = (xi • • • Xn) 1 / 71 (equal- weighting generating function). 

Under condition ( |60[ ) (in particular, in the latter examples) each process gjk(n(t)), 
t > can be rewritten as 

(62) g jk (M(Y(t))), t > 

for appropriate functions gjk on A™ < and all (j, k) G I 2 . Thus, in this case, the 



drift process fj in ([57]) can be written as 
fl(*)= ' 



2G(/i(t)) 
-1 



J2 D hi G(p(t))n h (t)iH{t) (<*w - m(t))(S a - w(*)) 



(h,i)ei 2 



(el 



2G( M (*)) 



(h,i,j,k)el 4 



tei 



-1 



2G( M (*)) 



XI *fc(M i ))/*(j)(*)A t (fc)(*)X] -M(2)(*))(<*w -McoW) 



= «(y(t)), f > o, 

where the function 2 : [0, oo)" -1 
1 



2G(M(y)) 



J is defined by 

^ ^(M( y ))M J (y)M fc (y)X (<5 j7 -M i (2 / ))(4/-M i (y)). 



U,k)eP- 



lei 



Finally, using the stationary distribution v of the spacings process described in the 
introduction, we define a normalized version u : [0, oo)™ -1 — > R of u by 



(63) 



u (y) : = u (y) - v i u ) = u(y) - 



[0,oo) r ' 



u{z) dv{z). 



Combining Theorem 1, ( 56 ) and the representation of the drift process g in terms 
of the spacings process Y, we can compare the value process V* of the portfolio 
generated by a permutation invariant function G G C to the value process of 
the market portfolio with the same initial value. Similarly, instead of (56), we can 
use 



log(^(i)/^(t)) = -log(G( / i(i))/G(^(0)))+ / (-flOO)da 

Jo 

and apply Theorem 1 with —u. All in all, we obtain the following corollary. 



Corollary 9. Suppose that G is permutation invariant, and that u in ( |63[ ) satisfies 
< Var u (u) = a 2 < oo, ||u||oo £ (0, oo) and 5 2 (u) > 0. Then for every initial 
spacing distribution k ofY(0) the ratios between the portfolio value V" (•) generated 
by G and the market portfolio value V""(-) satisfy the estimates 

P(V"(t)/V»(t) > C +(t)G( M (t))/G( M (0))) < 

P(V(t)/V"(t) > cr(t)G( M (0))/G( M (f))) < 
for all t,r,e > 0, where c^(t) :— exp[{r ± v(u)}t], j3 is given by Q and 

rax ^ 



dn 




r t i 


dv 


exp 

2 




dn 




r * i 


dv 


exp 

2 





C2 



,<5 2 (u) 



> 4e(e + OU/l + 



2£(e + cr 2 ) 2 ||u|| 



- 1 
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Example 1 . We now examine the statement of Corollary [9] in the cases of the port- 
folios generated by the diversity, the quadratic Gini coefficient, the Renyi entropy, 
the entropy and the equal-weighting generating function defined above. 

For the diversity-weighted portfolio, that is, the portfolio generated by the diver- 
sity, we compute 



(D jk G)( X ) = (i-p)(j2^) zf^r 1 + *j*(p - 1) ( £ * 



for all (j, fc) G I 2 . In addition, we note that the vector M(V(i)), which we will 
abbreviate by M(£) from now on, satisfies Mi(t) < M 2 (t) < ... < M n {t) for 
Lebesgue almost every t > with probability one. Indeed, this is a consequence of 
the fact that the set 

{t > : X^)(t) = for some 1 < i < j < n} 

is a Lebesgue zero set with probability one (see [55] for more details). Hence, 
we may conclude gjk(M(t)) — (DjkG)(M(t)) for Lebesgue almost every t > 
with probability one. Using the formula above for the partial derivatives Dj k G, 
(J, k) G I 2 and taking into account the latter observation, we obtain after some 
elementary computations: 



(EieiWW)") 

for Lebesgue almost every t > almost surely. Moreover, since by definition the 
market weights Mj(i), i G / are [0, l]-valued, we conclude using the Cauchy-Schwarz 
inequality that u(Y(t)) takes values in 0, ^tiW^El fo r Lebesgue almost every 
t > with probability one. Hence, the fluctuations of the relative performance of 
the diversity-weighted portfolio with respect to the market portfolio are controlled 
by the estimates of Corollary M with ||u||oo replaced by ai fo <5( u ) replaced 

by ("-Da-P) 

J n 

For the portfolio generated by the quadratic Gini coefficient a similar computation 
yields 

(D jk G)(x) = -S jk , (j,k)£l 2 

and, hence, 

- fvm , T^Mm? nZieAW)) 3 ) + (Ei £J (M(f)) 2 ) 2 

for Lebesgue almost every t > almost surely. Moreover, the inequality (Mj(i)) 2 — 
2(Mi(i)) 3 + {Mi{t)) A > for every t > and i G I, and the fact that the market 
weights are [0, l]-valued and sum up to one imply: 

< u(Y{t)) < 



2-il-n- 1 ) 2 

for Lebesgue almost every t > with probability one. Thus, the inequalities 
of Corollary M apply with ||ii||oo replaced by 2 -(i-7 t - 1 ) 2 an< ^ replaced by 

2 

l-i(l-n-!) 2 ' 
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In the case of the portfolio generated by the Renyi entropy we get 

o p— 1 p— 1 p— 2 

(D jk G)(x) = ■ > " - 6 jk 



P-l (E^?) 2 z2iei 



X ■ 



for all (j, k) E I 2 . An analogous computation to the case of the diversity- weighted 
portfolio yields here: 

P 



u(Y(t)) = 



21og(E ie /^(*) p ) 
ft PIP ^ Mi ^ P 2 ^ (Mi(r ' |VMfrf 



for Lebesgue almost every f > almost surely. Although the values of the process 
Mn(t), t > can be arbitrarily close to one, a simple analysis based on L'Hopital's 
rule shows that the values of the process t > are uniformly bounded 

for Lebesgue almost every t > with probability one, so that Corollary [9] applies 
in this case as well. 

In the case of the entropy-weighted portfolio, that is, the portfolio generated by 
the entropy, one computes 

{D ]k G){x) = -S jk —, (j,k) G I 2 , and 



u(Y(t)) = 



x 



3 

2 



for Lebesgue almost every i > almost surely. Moreover, the estimate 

|w(F(t))l ^ -2M n (t)logM n (t)' 

an analysis of the right-hand side of the latter inequality as M n (t) approaches one, 
and the inequality M n (t) > n^ 1 show that the values of the process |u(Y(f))|, 
t > are uniformly bounded for Lebesgue almost every t > with probability one. 
Hence, our Corollary [9] can be also applied in this case. 

Finally, for the egual- weighted portfolio, that is, the portfolio generated by the 
equal-weighting generating function, one easily checks 

(Dj k G)(x) = — - — (xi ■ ■ ■ x n ) 1 ^ n (- - 5jk\ , (j,k)el 2 , and 

, , . . n— 1 
u(Y(t)) = — 

for Lebesgue almost every t > almost surely. Thus, this is a trivial case and, 
although Corollary [9] applies, it does not give a meaningful estimate. 

4.2. Fluctuations of the market weights. On pages 46 and 52 of their survey 
on Stochastic Portfolio Theory [TS] , Fernholz and Karatzas pose the following open 
questions which we restate slightly for the models considered in this article. Con- 
sider an abstract rank-based equity market model with n companies as defined in 
the previous subsection and consider for any given time t > the ranked market 
weights: 

M(l)(*) < M(2)(*) < • < M(n)(*)- 

What can one say about the following objects: 
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(i) approximate laws of fim(t) and /Z( n )(i), 

(ii) fluctuations of the moving averages 

1 f T 

^ J H(k)(t)dt, k=l,...,n. 

(iii) In addition, the following question is of interest: What is the approximate 
deviation from 1 jn of the quantity 

— {amount of time the i-th market weight has rank k during [0, T]} . 

Our estimates on the rate of convergence to equilibrium will allow us to partially 
answer each of these questions. 

To answer question (ii), we recall from the previous subsection that at any given 
time each ranked market weight can be written as a time-independent continuous 
bounded function of the vector of spacings. Thus, our Theorem [l] can be applied 
directly to the moving averages in question (ii). It gives completely explicit esti- 
mates on the fluctuations of the latter, provided that one can compute the first 
two moments of the corresponding ranked market weight in equilibrium. This is 
in general a daunting task. In the following theorem we provide formulas for all 
moments of the ranked market weights in the Atlas model under their stationary 
distribution. Although not explicitly numerical, they can be effectively computed 
via a software such as Mathcmatica. This also gives a partial answer to question (i) 
for the Atlas model. Recall that the latter is the special case of the particle system 
in Q with Si = 5 > 0, 5 2 = . ■ ■ = S n = 0, <7i = . . . = a n = 1. 

Theorem 10. Consider the Atlas model with 8 being the drift of the lowest ranked 
particle. In equilibrium, the law of the ranked market weights is determined by the 
following Laplace transform: 



(64) t(0) := E 



exp 



where 

(i) <j) is the Laplace transform of e w with W being an Exponential random 
variable of parameter 25 /n, 

(ii) ipj is the conditional Laplace transform of 

-p \ n / 2S 



conditional on the value of f3, where ft is Beta(n — k + 1, k) distributed, V 
is uniformly distributed on (0,1), and ft, V are independent. 
In particular, in equilibrium, we obtain all moments of ixha by the formula 

i r°° 

E [(m) r ] = ^—YyJ o e r -\{e)de, r = i,2,.... 

Proof. We fix a k € {1, . . . , n} and recall the following result from Pal and Pitman 
[521 Theorem 8]. Let £i)£2i . . . , £n-i be independent Exponential random variables 
with respective parameters 

25 . 

— (n — i) , i = 1, 2, . . . , n — 1. 
n 
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Then take £o to be any random variable and set r\i := £o + £1 + . . . + 
i — 1, 2, . . . , n. Then, the following equality in law holds: 

M(fc) = — — , fc = 1,2, ...,n. 



Note that £o does not play any role, since it gets cancelled in the latter fraction. 
Thus, 



(65) - = ^ + l + ^ = 4 + l + S, 

M(fc) 



Hereby, 



fc — 1 fc— 1 n J — 1 



A fe = E exp (~E£0' Efe = ^ ex p(E^) 

j = l l=j j=k+l l=k 

are independent random variables. 

Next, let $i, i?2, • • • j denote i.i.d. Exponential random variables with some 
parameter a. Then the Renyi representation of the order statistics of i.i.d. Expo- 
nential random variables states that the random variables 

-#(»), * = 1,2,...,™- 1, 

are independent and Exponentially distributed with respective parameters a{m—i), 
i = 1, 2, . . . , m — 1. Hereby, < i?( 2 ) < . . . < #( m ) are the order statistics of the 
vector 1? 2 , . . . , #m)- We shall use this representation to express and in a 
symmetric way. 

Now, set m = n — k, a = 25 /n and define Aj := — i = 1,2, . . . ,m— 1, 

Ao := Then, we get the following equality in distribution: 

A,=£fc +i , 0<i<n-fc-l. 

Hence, it holds 

ro— fc n — fc n— fc i— 1 

E^ = E^E ex p(E A 0= Sfe 

i=l i=l i=l ;=0 

in distribution. Thus, is the sum of (n— fc) i.i.d. random variables. In particular, 

(66) i?( e -^) = (0 Q (0)r- fc , 

where (f) a is the Laplace transform of e 1 given by 

/>OG 

(67) 4> a {0)= / aexp ( - - 0e x ) dx. 

Jo 



The case of is a bit more convoluted. First, let T\,Ti, ■ . . ,T n be i.i.d. Expo- 
nential random variables with parameter a = 25 /n, which are independent of the 
?Vs. Setting Rj — cxp(-Tj) for j = 1, . . . , n, it is clear that each random variable 
Rj is distributed according to the Beta distribution Bcta(a, 1). Hence, we can write 
Rj = V X J a , j = 1, • • • , n with suitable i.i.d. uniformly on (0, 1) distributed random 
variables Ui, U2, ■ ■ ■ , U n . 
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Now, using the Renyi representation again we obtain the following identity in 
distribution: 



k— 1 —t,i\ fe — 1 d 1 / 7t \ l/ a n /rr \ !/ Q 

where T (1) < T (2) < . . . < T (n) , < JZ (2) < . . . < and C/ (1) < £7 (2) < . . . < 
[/(„) are the order statistics of the vectors (Ti, T 2 , . . . , T n ), (i?i, R 2 , . . . , R n ) and 
(Ui, U 2 , . . . , U n ), respectively. 

We now employ some known identities related to the Uniform distribution. First, 
we note that the vector 

\U{i),U(2) - J7(i), £7(3) - £7(2), . . • , £7(„) - £7( n _i), 1 - £/(„)) 

is distributed uniformly over the (n + l)-simplex {x e R n+1 : x, > 0, £. x i = !}j 
i.e. as Dirichlet(l, 1, . . . , 1). 

By the aggregation rule for the Dirichlet distribution the vector 

(£7(n-fe+i)) U( n - k+ 2) — £7( n _ fe+1 ), . . . , £7(„) — £7( n _i), 1 — £/(„)) 

has the Dirichlet(n — k + 1, 1, . . . , 1) distribution on the (fc + l)-simplex. 
Hence, by the usual Beta-Gamma algebra, we see that 

J} (£7(n-fc+2) — £7(n-fc+l)> • ■ ■ ) £7(n) ~ £7( n -l) , 1 ~ £7(n)) 

1 — U(n-k+l) 

is distributed as Dirichlet(l, 1, . . . , 1) over the 7c-simplex independently of U( n ^k+i)i 
which is distributed as Beta(n — k + 1, k). 

As a corollary, taking partial sums, we deduce that the law of the vector 

£7(n-fc+2) — £7(n-fc+l) £^(n-fc+3) — £^(n-fc+l) £7(n) — £7(n-fc+l) 



I - £7( n -fc+l) 1 — £7( n _fe+i) 1 — £7( n _fe_|_i) 

is the same as that of the order statistics of (k — 1) i.i.d. Uniform(0, 1) random 
variables V±, . . . , Vh~\ independent of /3 :— U( n -k+i)- 

Using the expression in ( 68 ) we obtain that Ak has the same law as 

fc-i . -s x i/c 

Hence, 

E (e- eA ») = E[^{6) k - 1 ] 
where ip-^ is the conditional Laplace transform of 



(1-0)^+0 



conditioned on /3. 

Hence, from ( [65] ) we get 

t(6) = E (e~ 9 /^) = e-°<l> a (0) n - k E[il^{0) k - 1 ], 
which leads to (64). 



28 



TOMOYUKI ICHIBA, SOUMIK PAL, AND MYKHAYLO SHKOLNIKOV 



To find the moments, we use the following fundamental identity: for any r > 0, 
we have 

9 r-l e -6/„ de = 

Replacing /j. by above and interchanging expectation and integral we get 

1 f°° 

This completes the proof. □ 




Unfortunately, the beautiful identities provided by the Atlas model do not extend 
to more general models. Asymptotic derivation of moments (when n tends to 
infinity) is possible in certain regimes due to an approximation by the atoms of the 
Poisson-Dirichlet distribution. Please see the article by Chatterjee and Pal [TT] for 
the details. Formulas for moments in the two-parameter Poisson-Dirichlet model 
can be found in Pitman and Yor [341 Proposition 17]. 

We now give an answer to question (iii) in the case of the particle system in ([2]) 
under the condition |3]). Recall that we assume = 1, i £ I throughout. 

Theorem 11. Let 

X{t) = (x 1 {t)-n- 1 Y J Mt),---,X n (t)-n- 1 Y J Mt)), t>{) 

be the centered version of the particle system in ^ and assume that ([3| holds. 
Then the process X is Markovian and possesses a unique invariant distribution v. 
Moreover, for every measure k which is absolutely continuous with respect to v and 
such that jj5 is square integrable with respect to v , one has for all t,r, e > the 
estimate 



P[ - I u(X(s))ds >r)< 



dv 



exp 

2 




1 max I — — ,4e(e + a 2 ) I J 1 -\ ' ^ 1 



Cp ' U 2 (m)' IV 2e(e + a 2 ) 2 ||u| 



for all bounded measurable functions u provided that the initial value X(0) is dis- 
tributed according to K, v(u) = and Varp(u) = a 2 . Hereby, Cp is a positive 
constant depending only on n and 5\, . . . ,S n (see ( |73[ ) for an explicit expression) . 
In particular, the latter estimate holds for functions of the form 

with a 2 — n~ 2 (n — 1), S(u) — 1 and = 1 — n^ 1 . 

Proof. 1) The Markov property and the existence and uniqueness of the invariant 
distribution of the process X were shown in Theorem 8 of |32j . Thus, we only 
need to prove the inequality in the statement of the theorem. To this end, we 
introduce for each vector x £ K™ a permutation ir(x) of the set {1, . . . , n} such that 
x ir(x)(i) ^ ' * " < x Tr(x)(n) holds. Since the process A is a diffusion process with state 
space 

fi = {i£l": Xx H \-x n = 0}, 
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it is a Feller process, the space C£°(.£j) is a core for its generator £ and on that 
space the generator is given by 

(69) (Cf)(x) = -^ Oij J ^x) + ^2fi <x) -i ({i -^-(x), xEft, 

i.j = l 1 J i=l 1 

where <Zjj = §ij — jli = Si — n z^jgi^j an( l / i s the composition of the 
projection of vectors in M. n onto f) and / (see chapter 18 of [25] for the details). 
Next, we define the cone 



f}< = {a; € fj : x± < • • • < £„} 



and introduce the mapping 9 : E" — > M™, which arranges the coordinates of a 
vector a; G R n in ascending order, as well as the mapping 

$: #< (K + ) n_1 , i4(i 2 -ii,...,i n -i»-i), 

which maps vectors in ,f)< to the corresponding vectors of spacings. We recall from 
Theorem 8 in [32] the following facts. The invariant distribution v of the process 
X is absolutely continuous with respect to the Lebesgue measure on Sj, its density 
is proportional to e~22fc=i *<=( e O)) Q <= anc j the process X is reversible with respect 
to v. 

2) In view of the results in step 1, as well as Theorem 3.1 in [T^ it suffices to show 
that the Poincare inequality 

(70) f f{x) 2 e- X«=l **( (*)W dx<C P f (-Cf)(x)f(x)e- **( e W) a » dx 

holds for a suitable constant Cp > 0, whereby the integration is performed with 
respect to the Lebesgue measure on the hyperplane fj. To this end, let £ be a 
positive definite symmetric n x n-matrix such that E 2 = (<Zy)i<j j< n . Then the 
same computation as on the top of page 64 in [3] , but with V replaced by SV, shows 



that inequality ( 70 1 is fulfilled provided that there exists a function V : 9j [1, oo) 



which belongs to H 1 ^; v) and satisfies 

CV 1 

( 71 ) -TT < 



V ~ 2C 



p 



almost everywhere on Sj. Hereby, H 1 (Sj; v) is the space of square integrable func- 
tions with respect to C>, whose gradient exists in the weak sense and is square 
integrable with respect to v. 

We claim that there is a < c < 1 such that the function 



(72) V(x) 



defined on has the desired properties. Indeed, V is a Lipschitz function and, thus, 
differentiable almost everywhere. Moreover, the condition < c < 1 shows that V 
belongs to i/ 1 (f);i>). In addition, it holds J2k=i ®k(®( x )) a k > by definition, so 
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that we have V > 1. Finally, we compute 



n 



2 

f=i 



2 n n 
i,i=l i=l 



^max \ II ~ ||2 



where we have set ao = a„ = and have written X max for the maximal eigenvalue 
of the matrix (a«)i<i *•<«.■ It is not hard to see that the eigenvalues of the latter 
are given by 1, ... ,1,0, so that X max = 1. Thus, for any c € (0, 1) the Poincare 



inequality (70) with the constant Cp — — ^an-j cTM P holds true. Taking the limit 



c t 1, we conclude that the Poincare inequality (70) is satisfied with 
(73) C ' ' 



p 



m 2 EieASi-n-'E^Sj) 2 ' 
This finishes the proof. □ 
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